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Abstract. A wavelet analysis has been performed on 80 intense gamma-ray bursts 
(GRBs) from the BATSE 3B catalog with durations longer than 2 seconds. The wavelet 
analysis applied novel features developed for edge detection in image processing and 
j> • this filtering process was used to extract a fit to the irregular GRB profile from the 

On ' background. A straightforward algorithm was subsequently used to identify statisti- 

cally significant peaks in this profile. The areas and FWHM of 270 peaks that were 
characterised as isolated were found to be consistent with lognormal distributions. The 
distribution of time intervals between peak maxima for all 963 identified peaks in the 
qq ' GRBs is also presented. 
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The recent significant discoveries of fading x-ray emission [1] have resulted in 
the first identification of radio and optical counterparts to GRBs [2,3]. The optical 
spectrum of one counterpart (GRB 970508) has yielded a value of z=0.835 [4]. 

GRBs produce an immense amount of energy (~ 10 51 ergs) which must be emitted 
by a medium with highly relativistic velocities (7 > 100). Most GRBs are highly 
variable with a variability scale much smaller than their overall durations and there 
is considerable debate as to the cause of the variability, which may be due to internal 
or external shocks [5-7]. 

GRBs generally have complex time profiles and previous work [8] has indicated 
that some parameters are consistent with lognormal distributions e.g. durations 
of GRBs and time intervals between peaks. BATSE bursts with multiple peaks 
were also analysed and the applicability of lognormal distributions in GRBs was 
confirmed [9]. Other techniques have also been used to profile the time structure 
of GRBs [10-12]. The results presented here on a large sample of BATSE GRBs 
were obtained using wavelets. 
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FIGURE 1. A wavelet transformation representation of GRB 920513 (BATSE trigger 1606) on 
three scales and the remaining low frequency signal. Note that the extreme values in the transform 
correspond to the edges in the signal and the signal can be reconstructed from just these extreme 
values [13]. This fact allows the removal of noise structure by removing the maxima and minima 
corresponding to the noise edges before the signal reconstruction [14]. 

DATA PREPARATION 



The data was obtained from the BATSE 3B catalogue. A subset of this catalogue 
was selected based on criteria used by Norris and co-workers [11], which gave a 
sample of 87 long bright bursts (T 90 >2 s and P256ms > 4.6 ph cur 2 s _1 ). The start 
and end times for each burst were identified by eye and a margin of 10 seconds 
added. Two sections of duration 30 s, one at -50 s from the start and the other 
at +20 s from the end were used to provide a linear interpolation for background 
subtraction. 

From the original sample of 87 bursts, a smaller sample of 80 background sub- 
tracted GRBs suitable for denoising (those with no data-gaps) were processed with 
a multiscale edge detection wavelet denoising routine [13-16]. 

Following on the wavelet denoising of the original GRB, a method was devised 
to identify the peak structures within the burst. Starting at the beginning of the 
burst, each peak was examined for minima on both sides which were separated from 
the maximum by more than a chosen significance level. If the search for minima 
failed on this peak it was rejected and the search continued. The algorithm was 
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FIGURE 2. The distribution of the areas of 270 peaks that were isolated at the 66% level (see 
text). The dashed line is a very acceptable lognormal fit to the data (see Table 1). 

designed so that the lowest minima were selected for each maxima. The analysis was 
extended to allow the identification of peaks which are well separated or isolated. 
The fraction of the total height which was above the higher minimum for every 
peak was calculated and then a threshold was applied on this fraction above which 
peaks were considered isolated. 

RESULTS 

Having identified peaks in the GRB time profiles, certain key characteristics of 
the data were examined. In particular the time interval between peaks (AT), the 
area under peaks (measured as total counts above background while the signal 
is above the height threshold for isolation) and the full width at half maximum 
(FWHM) were investigated. The distributions of each of these properties (Figures 
2 and 3) was tested for compatibility with the lognormal distribution. The results 
of the statistical tests are presented in Table 1. It should be noted that the peak 
area and FWHM are quite compatible with the lognormal distribution. 

CONCLUSIONS 

A comprehensive profile of the temporal properties of intense GRBs has now 
been provided by the wavelet method. Many properties of GRBs are consistent 
with lognormal distributions: 
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FIGURE 3. The FWHM distribution for the same peaks as in Figure 2. Again the dashed line 
is a lognormal fit to the data. 

• The bimodal duration distribution of GRBs is consistent with two lognormal 
distributions [8]. 

• The area of the peaks is consistent with a lognormal distribution. 

• The FWHM of the peaks is consistent with a lognormal distribution. 

• The time intervals between peaks in GRBs as measured with this technique 
does not seem to be consistent with a lognormal distribution but important 
selection effects have not been considered. 

These results are in good agreement with previous work [8,9]. These results pro- 
vide an alternative description of GRB profiles to that of the peak-aligned stretched 
exponentials [10-12]. The data presented here also provide considerable constraints 
on models that attribute the variability to interactions with an external medium 
or to an internal origin that reflects the activity of the central engine. 



TABLE 1. Summary of goodness of fit tests for lognormal 
distributions, x\a ^ the 2a confidence limit of the reduced 
X 2 value (x 2 jv) with v degrees of freedom. 
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FIGURE 4. The distribution of time intervals between the 963 peaks identified in 80 GRBs. The 
dashed line indicates a candidate lognormal fit. The departure of the fit at short time intervals 
is strongly influenced by the limited resolution of the data [17]. 
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